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A detailed microscopic and quantitative description of the electronic and magnetic properties of 
Gd 3+ -doped YC 02 Z 1120 single crystals (Yi_ x Gda;Co 2 Zn 2 o: (0.002 < x < 1.00) is reported through 
a combination of temperature-dependent electron spin resonance (ESR), heat capacity and dc mag¬ 
netic susceptibility experiments, plus first-principles density functional theory (DFT) calculations. 
The ESR results indicate that this system features an exchange bottleneck scenario wherein various 
channels for the spin-lattice relaxation mechanism of the Gd 3+ ions can be identified via exchange 
interactions with different types of conduction electrons at the Fermi level. Quantitative support 
from the other techniques allow to extract the exchange interaction parameters between the local¬ 
ized magnetic moments of the Gd 3+ ions and the different types of conduction electrons present at 
the Fermi level (J/ s , Jj p and Jfd)- Despite the complexity of the crystal structure, our combina¬ 
tion of experimental and electronic structure data establish GdCo 2 Zn 2 o as a model RKKY system 
by predicting a Curie-Weiss temperature 6c = —1.2(2) K directly from microscopic parameters, 
in very good agreement with the bulk value from magnetization data. The successful microscopic 
understanding of the electronic structure and behavior for the two end compounds YCo 2 Zn 2 o and 
GdCo 2 Zii 2 o means they can be used as references to help describe the more complex electronic 
properties of related materials. 

PACS numbers: 76.30.-v, 71.20.-b, 76.60.Es 


I. INTRODUCTION 


The understanding of the physical properties of com¬ 
plex materials with cage-like structures, such as the fam¬ 
ily RT 2 Zn 2 o (R = rare earth, T = transition metal), have 
attracted the attention of many researchers focused in 
condensed matter physics. Among other aspects, these 
intermetallic cage compounds generate interest due to 
the different types of electronic and magnetic b ehav iors 
governed by “naturally diluted” rare earth ionsP^l The 
magnetic versatility associated with the 4/ electrons goes 
from weakly correlated Pauli-like paramagnetic behav¬ 
ior (Lu 3+ ), to hybridization with conduction electrons 
(Yb 3+ ) and peculiar interactions of local magnetic mo¬ 
ments (Gd 3+ ) with conduction electrons, as a few exam¬ 
ples. 

Despite the apparent complexity (184 atoms per con¬ 
ventional unit cell arranged in different types of cages), 
their crystallographic structure can be broken down into 
surprisingly simple sub-units, which allows clean analyzes 
and interpretations of experimental results, and sets this 
family as an excellent model system for several physical 
problems. They adopt a cubic CeCr 2 Al 2 o-type structure 
(space group: FdSm )P in which the R and T ions occupy 
their own unique crystallographic sites (8a and 16d, re¬ 
spectively). The Zn ions form the cage structure by occu¬ 
pying three inequivalent crystallographic sites (96<?, 48/ 
and 16c). If we consider the sub-structure in terms of the 
nearest neighbors and the next nearest neighbors, the R 
and T ions are fully surrounded by shells formed by Zn 
ions, leaving a shortest R-R spacing of r ss 6 A. The R 


ions are thus isolated in Frank-Kasper cages formed by 
16 Zn ions as exemplified in Fig. [I] 

a) 
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FIG. 1: (a) Conventional unit cell of YCo 2 Zn 2 o (space group 
Fd3m). Y, Co and Zn ions are represented by red, blue and 
gray balls, respectively. The CN-16 Frank-Kasper polyhe¬ 
dron around the Y ions and the polyhedron around Co ions 
are highlighted in pink and violet, respectively, (b) Typical 
crystal of YCo 2 Zn 2 o showing a [111] facet. 

The observation of radically different magnetic behav¬ 
iors such as the contrast between low-temperature anti¬ 
ferromagnetic order (T)v ~ 5.7 in GdCo 2 Zn 2 o and 
the high temperature ferromagnetic order {Tq ~ 86 K]P 
in GdFe 2 Zn 20 hints at peculiar magnetic couplings gov¬ 
erned by the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction. In the case of Y-based compounds, largely 
different types of behavior are also found. YFe 2 Zn 20 has 
been described as a “nearly ferromagnetic Fermi sea” 
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compound 5 because it is near the Stoner limit, in contrast 
to YCo 2 Zn 2 o with more conventional metallic behavior. 
Moreover, within this family all the Yb-based compounds 
described so far have presented heavy fermion behavior 
as evidenced by the enhanced Sommerfeld coefficients, 
reaching 7 ss 7900 mJ/mol.K 2 for YbCo 2 Zn 2 o- All of 
these observations point to the need to investigate the 
electronic structure,® particularly around the Fermi level, 
and describe the electronic interactions in detail to better 
understand the electronic and magnetic behaviors. 

In order to conduct this task in a tractable manner, 
we have chosen an initial focus on the weakly correlated 
compound YCo 2 Zn 2 o that features Pauli like paramag¬ 
netism, metallic transport and a Sommerfeld coefficient 
of 7 = 18.3 mJ/mol.K 2 pEl as an appropriate host for a 
microscopic study using Gd 3+ ions as an electron spin 
resonance (ESR) probe. ESR of rare earth ions diluted 
in metallic hosts is a useful local technique to investi¬ 
gate microscopic properties of materials, since it directly 
probes the localized magnetic mome nts and the nature of 
the interactions with their neighbors !- 8 9 The metallic and 
nonmagnetic YCo 2 Zn 2 o host doped with Gd is an excel¬ 
lent model system to study the Gd 3+ spin-lattice relax¬ 
ation, associated with the conduction electrons (ce) spin- 
flip scattering mechanism due to the exchange interaction 
between the localized ma gnetic moment and the ce. The 
Hasegawa-Korringa modef^ 11 for the spin-lattice relax¬ 
ation has been carefully discussed and applied in previous 
studies of ESR for Gd 3+ in the intermetallic compounds 
like LaAlJ 12 } LuAl ^ 3 and in elemental AlP® We have 
recently applied the same technique to investigate the 
electronic structure of the superconductor YIn 3 P® 

In this work we show that the ESR spectra of 
Gd 3+ in YCo 2 Zn 2 o ( 0.001 < x < 1 . 00 ) presents a 
Gd 3+ concentration-dependent thermal broadening of 
the linewidth and 3 -shift. This reveals the existence of 
the exchange bottleneck effect in this compound, that can 
be tuned by the concentration of Gd 3+ . By combining 
the ESR results with heat capacity, magnetic suscepti¬ 
bility and band structure calculations, we extract the 
exchange parameters of the interaction between Gd 3+ 
and the s, p, and d ce present at the Fermi level of 
YCo 2 Zn 2 o- We are then able to establish a clear cor¬ 
relation of these microscopic parameters with the RKKY 
interaction. This in turn offers a better understanding 
of the peculiar “Fermi sea” present in the system, which 
has lead to the magnetic anomalies found in this family. 


II. EXPERIMENTAL AND COMPUTATIONAL 
DETAILS 

Batches of Y 1 _ x Gd x Co 2 Zn 2 o (0.001 < x < 1.00) single 
crystals were grown by the self-flux method^® 0 using ex¬ 
cess Zn. The starting reagents were 99.9% Y, 99.9% Co, 
99.9% Gd and 99.9999% Zn (Alfa-Aesar). Initial ratios 
of elements were 1:2:47 for the pure ternaries Y:Co:Zn 
and Gd:Co:Zn, or l-x:x:2:47 for the pseudo-quaternaries 


Y:Gd:Co:Zn, based on previously reported growths of the 
ternary compounds.®The elemental mixtures were sealed 
in an evacuated quartz ampoule and placed in a box fur¬ 
nace for the temperature ramping. Crystals were grown 
by slowly cooling the melt between 1100 °C and 600 °C 
over 100 h. At 600 °C the ampoules were removed from 
the furnace, inverted and placed in a centrifuge to spin 
off the excess flux. The separated crystals are typically 
polyhedral, ~ 3 mm or larger and manifest clear, trian¬ 
gular [111] facets (Fig. [l]). The Gd concentrations were 
estimated based on the effective moments per formula 
unit extracted from fits of magnetic susceptibility mea¬ 
surements. Powder x-ray diffraction on crushed crystals 
was used to ascertain the CeCr 2 Al 2 o-type structure^ as 
exemplified for YCo 2 Zn 2 o in Fig. [ 2 ] The refined lattice 
parameter of a = 14.042(4) A is in good agreement with 
the literature.® The inset shows that the refined lattice 
parameter increases linearly with Gd concentration, as 
expected by Vegard’s law.®** 



FIG. 2: XRD pattern of YCo 2 Zn 2 o- Inset, Lattice parameter 
for the Yi_ a! Gd iI ;Co 2 Zn 2 o ( 0.001 < x < 1 . 00 ) compounds. 

For the ESR experiments, single crystals were crushed 
into fine powders of particle size greater than 100 /iin, 
corresponding to average grain size ( d ) larger than the 
skin depth (<5), A = d/6 > 10. We have noticed that ex¬ 
periments on as-grown single crystals should be carried 
out with caution, because strong resonances similar to 
those reported by Ivanshin et al are frequently observ¬ 
able, and likely due to residual Co surface contamination 
since that particular signal disappears after removing the 
as-grown crystal surfaces. The X-Band (v k, 9.4 GHz) 
ESR experiments were carried out in a conventional CW 
Bruker-ELEXSYS 500 ESR spectrometer using a TE 102 
cavity. The sample temperature was changed using a 
helium gas-flux coupled to an Oxford temperature con¬ 
troller. The specific heat (C p ) and magnetic susceptibil¬ 
ity (x = M/H) measurements were performed on Quan¬ 
tum Design PPMS and SQUID-VSM platforms, respec¬ 
tively, using their standard procedures. 
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The ground state crystal structures were calculated 
using spin-polarized first-principles density functional 
theory (DFT), using the PBEsol exchange-correlation 
functional.^ The Kohn-Sham equations were solved us¬ 
ing the projector augmented plane-wave (PAW) method 
as implemented in the VASP code. 2122; The PAW 
atomic reference configurations are: 4s 2 4p 6 4d 1 5s 2 for 
Y, 5s 2 5p 6 4/ 7 5d 1 6s 2 for Gd, 3s 2 3p 6 3<T4s 2 for Co and 
3p 6 3d 10 4s 2 for Zn, where only electrons treated as valence 
electrons are explicitly enumerated. The energy cut-off 
in the plane-waves expansion is 507.5 eV, where the total 
energy has been converged tol meV/unit cell. All struc¬ 
tural parameters, lattice constants and atomic positions 
for each calculated compound have been optimized by si¬ 
multaneously minimizing all atomic forces and stress ten¬ 
sor components via a conjugate gradient method. Suc¬ 
cessive full-cell optimizations adapting basis vectors have 
been conducted until the unit cell energies and structural 
parameters were fully converged. Brillouin-zone integra¬ 
tion has been performed on a Monkhorst-Pack 12 x 12 x 12 
k-point grid with a Gaussian broadening of 0.01 eV for 
full relaxation (ionic forces are converged to 0.1 meV/A). 
Then, we used the relaxed crystal structures to calculate 
the total and partial density of states (DOS), the dis¬ 
persion relations and the Fermi surfaces using the full- 
potential augmented-plane wave method with local or¬ 
bitals!^! The muffin-tin (MT) radii of Y, Gd, Co and 
Zn are set to Amt = 2.6,2.8138,2.2804 and 2.2804 a.u., 
respectively. The parameter Amt | G + k| max governing 
the number of plane waves has been converged to 9.0. 
The irreducible wedge of Brillouin zone is sampled with 
a 20 x 20 x 20 uniformly spaced k-point grid. Fermi sur¬ 
faces were plotted using the XCrysDen package.^ 


III. EXPERIMENTAL RESULTS 


Figure [s] displays the Gd 3+ ESR spectra in 
Yi_ x Gd a; Co 2 Zn 2 o for x ~ 0.002 and x « 0.21 at 4.4 K 
and 12.5 K, respectively, and microwave power of P ^ ~ 
2 mW. These ESR spectra show different resonance mag¬ 
netic fields than that of Gd 3+ in insulators, for which 
a resonance field of H 0 = 3386(4) Oe and g -value of 
1.993(2) are well established.^ It is evident that for the 
low concentration sample (x = 0.002) the resonance is 
shifted toward a lower field (higher ^-values) compared 
to that of the higher concentration sample (x = 0.21). 

The observed ESR spectra of Gd 3+ localized magnetic 
moments in Yi_ a; Gd a ;Co 2 Zn 2 o will be analyzed according 
to the generally accepted approach where, at resonance, 
the microwave absorption in a metal is given b y the 
Dyson theory in the diffusionless limit, A/B « 2.6P^ 27 
In this limit, for particles larger than the skin depth, the 
ESR spectra reduce to a simple admixture of absorption 
(x") and dispersion (%') of Lorentzian lineshapesP^ 27 ! 
The derivative of this admixture is given by: 



H(Oe) 

FIG. 3: Gd 3 + ESR spectra of Y 1 _ a .Gd !C Co 2 Zn 2 o for x « 0.002 
at T = 4.4 K and x ss 0.21 at T = 12.5 K for a microwave 
power of Ppu 2 mW. 
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* = (H 0 -H) le T 2 , 

where H 0 and H are the resonance and the applied 
fields respectively, q e is the electron gyromagnetic ra¬ 
tio, T 2 the spin-spin relaxation time, a the admixture 
of absorption (a = 0) and dispersion (a = 1) and \o 
the paramagnetic contribution from the static suscepti¬ 
bility. It is usually accepted 2 ^ that for diluted magnetic 
moments in a metallic host T \ ps T 2 , where Tj is the spin- 
lattice relaxation time. 2 ® 2 ® Therefore, the fitting of the 
experimentally observed ESR absorption lines to Eq. 1 
allows the extraction of the two most relevant ESR pa¬ 
rameters, i.e., the g -value from the resonance condition, 
hv = gg B H 0 , and the linewidth AH = l/q e T 2 . 

Figure |4] displays the T-dependence of the Gd 3+ ESR 
linewidth, AH, in Yi_ a; Gd. r Co 2 Zn 2 o for 0.002 < x < 
1.00 at a microwave power of « 2 mW. The broad¬ 
ening of AH at low temperature for the high concentra¬ 
tion samples is presumably originated by the interaction 
between randomly distributed Gd 3+ magnetic moments, 
that cause an inhomogeneous local field. However, for 
the stoichiometric GdCo 2 Zn 2 o this disorder should be ab¬ 
sent and the low temperature broadening of AH less pro¬ 
nounced (black symbols in Fig. [4]). This allows the drop 
in the magnetic susceptibility at T/v to become evident in 
the ESR intensity (see Fig. [7j) . The high T-dependence of 
AH follows the linear behavior AH = a + bT where the 
a parameter represents the residual linewidth, AHq, and 
b = d(AH)/dT the Korringa-like relaxation rate. The 
extracted a and b values are given in Table I, together 
with the obtained p-values at T ss 10 K for the studied 
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FIG. 4: T-dependence of the Gd 3+ ESR linewidth, AH, in 
Yi_ a; Gd :l ;Co2Zn2o for 0.001 < x < 1.00. 

Figure [5] displays the Gd concentration dependence of 
the g-shift (A g = g - 1.993(2)) and in the inset the ther¬ 
mal broadening of the linewidth, b. The general trends 
of the data presented in Fig. [5] is characteristic of an ex¬ 
change bottleneck phenomenon, where the ce relaxation 
to the Gd 3+ localized magnetic moment (Overhauser re¬ 
laxation) overcomes the ce spin-lattice relaxation. 



x 


FIG. 5: Gd concentration dependence of the g-shift (A g = 
g - 1.993(2)) and in the Inset the thermal broadening of the 
linewidth, b, for the Yi_ x Gda;Co2Zn2o system. The dashed 
lines are guides for the eye. 


7 = 18(3) mJ/mol.K 2 and a Debye temperature of 0£> 
= 370(7) K for the YCo2Zn 2 o compound. Our obtained 
value of 7 is in agreement with the previously reported 
value of 18.3 mJ/mol.K 2 .® The inset of Fig. ml zooms in 
on the weak T-dependence of the magnetic susceptibility, 
X(T), for YCo2Zn 2 o- The dome-like feature at T « 120 K 
may be due to T-dependent d-ce at the Fermi surface (see 
Fig. 10 below) and the small upward tail at low tempera¬ 
ture (which has no influence in the analyzes that follow) 
due to residual rare-earth magnetic impurities from the 
99.9% Y reagents used in our samples. 



FIG. 6: Low temperature specific heat, C p /T, showing 
a linear behavior with a Sommerfeld coefficient of 7 = 
18(3) mj/mol K 2 and a Debye temperature of 0u = 370(7) K 
for YC02Z1120. The inset zooms in on the weak T-dependence 
of the magnetic susceptibility, x(T), for YCo 2 Zn2o 


Fig. [7| presents the T-dependence of the magnetic sus¬ 
ceptibility of GdCo 2 Zn 2 o where the antiferromagnetic or¬ 
der is seen at a Neel temperature of T/v ~ 5.7 K, with an 
effective magnetic moment g e ff = 8.1(2) g B , comparable 
to that of Gd 3+ ions = 7.94 g B ) and a paramag¬ 

netic Curie temperature 6q = —0.7(4) K. Note that the 
paramagnetic Curie temperature 9c is the same as Tq in 
the Weiss molecular field theory. Fig. [7] also shows that 
the ESR signal integrated intensity follows the trends of 
the magnetic susceptibility, crossing the paramagnetic- 
antiferromagnetic transition at approximately the same 
Neel temperature. This is in itself a rare observation 
because usually, around the ordering temperature of a 
magnetic transition, the ESR signal is lost due to a strong 
broadening of the resonance. 


We now focus on bulk thermodynamic measurements 
which, together with the band structure calculations pre¬ 
sented in the following section, provide support for a 
proper quantitative analysis of the ESR results. Figure[6] 
shows the low temperature linear behavior of C p /T as 
a function of T 2 leading to a Sommerfeld coefficient of 


The magnetic susceptibility for all our 
Yi_ a ,Gd a; Co2Zn2o samples were fitted to a Curie- 
Weiss law using the effective magnetic moment of 
g e ff = 7.94 g B for the Gd 3+ ions. From these fittings 
the Gd concentrations were estimated and their values 
are listed in Table |l| 
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FIG. 7: T-dependence of the dc magnetic susceptibility and 
the Gd 3+ ESR intensity for GdCo 2 Zn 2 o. Both experiments 
show the magnetic paramagnetic-antiferromagnetic transition 
at Tn ss 5.7 K. The inset shows the inverse magnetic suscep¬ 
tibility and linear fit results. 


IV. ANALYSIS AND DISCUSSION 

The low-T linear behavior of C p /T = 7 + 0T 2 for 
YCo2Zn 2 o of Figure [6] leads to a Sommerfeld coeffi¬ 
cient 7 = 18(3) mJ/mol K 2 and a Debye temperature 
Qd = 370(7) K. In the Fermi liquid model the Som¬ 
merfeld coefficient is given by 7 = (2/3)7 T 2 k^r]F, where 
r/p is the total density of states (DOS) per formula unit 
(f.u.), spin and eV at the Fermi level. Thus, we estimate 
t]f = 3.8(8) states/f.u., spin and eV for YCo 2 Zn 20 . 

A Pauli-like paramagnetic susceptibility for YCo 2 Zn 2 o 
can be estimated at high-T (310 K). The data in the inset 
of Figure [6] presents, after correction by the core diamag¬ 
netism of YCo 2 Zn 2 o ( Xdia = —2.3 x 10 -4 emu/mol), a 
lower limit magnetic susceptibility value at T sa 310 K 
of xp = Xo — Xdia = 0.312 x 10 -3 emu/mol which is 
slightly smaller than the value reported by Jia et al.® 
Notice that the diamagnetism of the cage structure has 
not been considered. 

Once again within the Fermi liquid model, the Pauli- 
like paramagnetic susceptibility is given by xp = 

Then, using the experimental Pauli paramagnetic suscep¬ 
tibility for YCo 2 Zn 2 o we estimate rjF = 4.8(9) states/f.u., 
spin and eV as a lower limit for the DOS. However, 
one has to consider the possibility of an exchange- 
enhanced magnetic susceptibility, i.e., x = Xo/(l ~ £), 
where / accounts for the electron-electron exchange en¬ 
hancement. Nevertheless, the estimated value of rjp = 
4.8(9) states/f.u., spin and eV is, within the accuracy 
or our experiments, comparable to the value obtained 
from the Sommerfeld coefficient. Thus, in our analysis 
we shall ignore the electron-electron exchange enhance¬ 
ment for YCo 2 Zn 2 o. 

From the high-T data of Figure [7] (100 K - 300 K) 
for GdCo 2 Zn 2 o we obtain a small negative value for the 


TABLE I: Gd concentrations, (/-values, residual linewidths, 
a, and thermal broadening of the linewidths, b , for the 
Yi_ a; Gda;Co 2 Zn 2 o system. 


Cone. 

(/-value 

Measured b 

Calculated b 

X 

(10 K) 

(Oe/K) 

(Oe/K) 

0.002 

2 . 012 ( 1 ) 

8 . 6 ( 6 ) 

8.4(3) 

0.005 

2.006(1) 

7.1(6) 

4.3(3) 

0.085 

1.994(1) 

7.2(6) 

0.02(3) 

0.21 

1.992(1) 

6 . 1 ( 6 ) 

0.02(3) 

0.56 

1.991(2) 

4.8(6) 

0.08(5) 

1.00 

1.991(2) 

1.9(6) 

0.09(5) 


Curie-Weiss temperature, dc = -0.7(4) K, as expected 
for this antiferromagnetic material. This is in contrast 
to the small positive value reported in previous studies.® 
In order to complement the experimental data analy¬ 
sis and provide details about the electronic structure, we 
have performed first-principles band structure calcula¬ 
tions for the pure compounds GdCo 2 Zn 2 o and YCo 2 Zn 2 o- 
The calculated lattice constant of the fully relaxed crystal 
structures are 13.7738 A and 13.7564 A for GdCo 2 Zn 2 o 
and YCo 2 Zn 2 o, respectively, obtaining an absolute rela¬ 
tive error of 2%, as it is expected from local and semi¬ 
local functionals. For GdCo 2 Zn 20 we have obtained 
that the magnetic stabilization energy Aafm = Ea fm — 
Efm = —2 meV/unit cell, therefore the ground state is 
antiferromagnetic, with a calculated local magnetization 
of 7.03/rB/Gd ion and zero local magnetization for the 
Co ions. On the other hand, YCo 2 Zn 2 o converged to a 
non magnetic ground state. A previous band structure 
calculation^ with local and semi-local functionals with¬ 
out including explicitly the Gd /-states into the valence 
window and with a tiny k-point grid, found for the un¬ 
relaxed crystal structure of GdCo 2 Zn 20 a slightly larger 
local magnetic moment of 7.25/is/Gd ion. However, it 
is worth noting that the calculation of the magnetic mo¬ 
ment is very sensitive to the number of k-points and those 
values were not fully converged.® 

Figure [8] shows the calculated dispersion relations for 
these two systems. For GdCo 2 Zn 2 o the valence bands are 
built up mainly from Zn d-states and a small contribution 
of Zn s- and p-states between -11 and -6.5 eV for both 
systems. The remaining valence bands result from the 
hybridization of Gd d- and /-states (the latter localized 
between -3.3 and -2.8 eV), Co d-states and Zn s-, p- and 
d-states. The conduction bands are also built up from 
the hybridization of Gd d- and /-states, Co d-states and 
Zn s-, p- and d-states, with the Gd /-states localized 
between 0.3 and 1.0 eV. On the other hand, the bands 
of YCo 2 Zn 2 o are built up similar to the previous system 
with the obvious absence of /-states, in this case Y d- 
states contribute mostly to the upper valence bands and 
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FIG. 8: Calculated dispersion relations for the GdCo2Zii2o 
(upper) and YC02Z1120 (lower) systems. Highlighted in colors 
are the four conduction bands crossing the Fermi level; band 
1 (yellow), band 2 (green), band 3 (cyan) and band 4 (blue). 
The eigenvalues are shifted with respect to the Femi level, 
which is indicated by a red line. 


the conduction bands. As can be observed in Fig. [HJ 
the topology of the bands is almost identical for the two 
systems, especially the four conduction bands that cross 
the Fermi level. The similarity of these four bands is more 
lively appreciable in the branches of the Fermi surface, as 
seen in Fig. [9} The first branch is formed by eight pockets 
along the Q direction of the first Brillouin zone. These 
pockets are formed from the contributions of Gd(Y) re¬ 
states and Co d-states. The second branch of the Fermi 
surface has eight connected structures along the A, V, Q 
and A directions that resemble a six-arm starfish. These 
starfish are built up from Gd(Y) d-states, Co d-states and 
Zn p and d-states. The third branch has one sphere at the 
center of the first Brillouin zone, which is made up from 
Zn s-states, and six structures that have the appearance 
of mushrooms with the stem along the A direction. These 
six mushrooms are formed from Gd(Y) d-states and Co d- 


states. The fourth branch is constituted of six lenses with 
their surfaces perpendicular to the A direction. These 
lenses are built up from Co d-states and Zn p-states. 



FIG. 9: Calculated Fermi surface for the GdCo2Zn2o and 
YCo2Zn2o systems. The four branches for GdCo2Zn2o cor¬ 
responding to band 1, band 2, band 3 and band 4 are shown 
in (a), (b), (c) and (d), respectively. In the same manner, the 
four branches of YCo2Zn2o are shown in (e), (f), (g) and (h). 


For GdCo2Zn 2 o our calculations estimate a total DOS 
at the Fermi level of 3.03(1) states/f.u., spin and eV 
(Fig. 10). Similarly, a total DOS at the Fermi level 
of 3.02(l) states/f.u., spin and eV was estimated for 
YCo 2 Zn 2 o (Fig. 10), which is comparable to the above 
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values obtained experimentally. Also, our calculated 
values are in relatively good agreement with a previ¬ 
ous first-principles study with a much less dense k-point 
grid.® These results support the previous statement that 
YCo 2 Zn 2 o may be considered as an intermetallic com¬ 
pound with negligible electron-electron correlations, i.e, 

X = Xo/(l - 0 with £ << 1. 

With all these details in mind, we can now return 
to the ESR analysis. The exchange interaction, H = 
-J fs Sf . s ce between the localized 4/-electron spin of 
Gd 3+ , Sf, and the ce of the YCo 2 Zn 2 o, s ce , yields an ESR 
5 -shift, Agj^l and thermal broadening of the linewidth, 
b, (Korringa rate) 1 !] given by: 


A g = Jf s {0)r] F 


( 2 ) 


and 


d(AH) 

dT 


7T k B 
g^B 


Jf S ( °)v 2 f = 


7T k B 
9dB 


(A g? 


(3) 


where Jf s ( 0) is the effective exchange parameter in 
the absence of ce momentum transfer, i.e., ( Jf s (q)) F = 
,// s ( 0 )P® fj F i s the “bare” density of states for one spin 
direction at the Fermi surface, k F is the Boltzmann con¬ 
stant, p F is the Bohr magneton, and g is the Gd 3+ g- 
value. 

The above equations are normally used in the analysis 
of the ESR data in the limit of very diluted rare earths 
and concentration-independent g and b parameters, i.e., 
in the non-bottleneck regime, absence of ce momentum 
transfer and single band compounds.®'Using the data in 
Table I and Eq. 3 one can observe that the correlation 
between 5 -shift, A 5 , and b is only verified for the lowest 
Gd concentration samples. Moreover, Fig. 0 Fig. 0 and 
Table [I] show clear concentration-dependent ESR param¬ 
eters, so our data must be analyzed in a different manner. 

The change of the Gd 3+ 5 -shift from positive values 
(at low Gd concentrations) to negative ones (at high Gd 
concentrations) and the x-dependence of the Gd 3+ ESR 
thermal broadening of the linewidth, 6 , (see Fig. 0 and 
Fig. 0) lead us to conclude that the relaxation of the 
Gd 3 “^ions to the lattice is processed via an exchange 
interaction, J(5/.s ce ), between the Gd 3+ localized mag¬ 
netic moment and different types of ce at the Fermi level. 

In a multiband approximation the 5 -shift, A 5 , and 
thermal broadening of the linewidth, d, are given by: 


A 5 = Ag fs + Ag fp + Ag fd 


= Jf s (0)ri Fs - Jf P (0)i] Fp + Jfd(0)VF d 


(4) 
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FIG. 10: Calculated total and partial DOS for the 

GdCo 2 Zii 2 o and YCo 2 Zn 2 o systems. The Fermi level is in¬ 
dicated by a dotted line. 


b = ^[F s Ag%+F p Ag% + F d Ag%] 

gg B 

(5) 

= ~^^[Fa{Jf a (Q)) Fr )F. + FpJfp(fy r l Fp + FdJfd(0) r l Fd ]i 


where k F is the Boltzmann constant, ps the Bohr 
magneton, and 5 the Gd 3+ 5 -value; Jfi(0) (i = s,p,d) 
are the effective 5 = 0 components of the exchange in¬ 
teraction between the Gd 3+ 4/ magnetic moment and 
the s-, p- and d-type ce; rj Fi (i = s,p,d) the partial 
bare DOS (states/f.u., spin and eV) at the Fermi level 
of the s-, p- and d-type ce; (Jf s (q ))_f is the average over 
the Fermi surface of the square of the 5 -dependent effec¬ 
tive exchange parameter in the presence of ce momen¬ 
tum transfer, 5 = \k out - k in |, i.e., ( Jf s (q)) F ^ J/ S ( 0)®1 
F s = 1, F p = 1/3 and F d = 1/5 are factors associated 
with the orbital degeneracy of the unsplit (no crystal 
field effects) bands at the Fermi level, respectively. The 
5 -dependence of the exchange interaction with the p-type 
and d-type ce will be considered constant over the Fermi 
surface, i.e., (J f P , d {q )}f = Jfp,d( 0 ) (see below). 

In Eq. 4 we have considered that the contribution to 
A 5 due to the exchange interaction with s- and d-type 
ce are positive (atomic-like) and that with p-type ce is 
negative (covalent-like) P® 

Due to the strong spin-orbit coupling of p- and d-type 
ce compared to that of the s-type ce, we assume that only 
the s-type ce are capable of experiencing the bottleneck 
effect. Hence, we can consider that the contribution to 
the ESR parameters, 5 -shift and b, of the s-type ce are 
negligible in the highly concentrated samples. With these 


and 























assumptions, for x = 1 (GdCo 2 Zn 2 o, extreme bottleneck ) 
Eq. 4 reduces to 


Ag = -0.002(2) = -Jf p {0)r)F p + Jfd(0)r]F d , ( 6 ) 
and Eq. 5 reduces to (see Fig. [l0]b) : 


b = 1.9(6) Oe/K = 0)vl p + F d J 2 fd ( 0)r£j.(7) 
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In general, the g-dependence of the exchange param¬ 
eters Jf d (* = s,p,d) cannot be disregarded, but in 
the case in which the relaxation rate, b, scales or is 
slightly larger than the expected value from the g-shift 
(b = ( 7 rfcs/g/rs)(Ag) 2 ) we can neglect the g-dependence 
of the Jf p and Jf d exchange parameters. In our case, we 
have that 1.9 Oe/K > (~ 2.34 x 10 4 Oe/IC)x(-0.002 ) 2 « 
0.1 Oe/K in agreement with Table I. 

From Fig. 10 for GdCo 2 Zn 2 o, we have q Fd = 1.09(1) 
states/f.u., spin and eV and r) Fp = 0.59(1) states/f.u., 
spin and eV. Then, using Eqs. 6 and 7 we estimate 
J fd ( 0) = 10(5) meV and J fp ( 0) = 22(6) meV. 

Conversely, in the not bottlenecked regime (lowest 
Gd concentration, x w 0.002), and from Fig. [lO] for 
YCo 2 Zn 2 o, Eqs. 4 and 5 reduce to: 


material with a complex Fermi surface is thus expected 
to be an almost impossible task. The equation that rep¬ 
resents a generalized form of the RKKY interaction is 
given by23 


Jrkky 



/k - /k+q 

e(k + q)-e(k)’ 


( 10 ) 


where /k = 0 (fc F —|k|) and /k+ q = 0 (fc F — |k + q|) are 
the step functions that come from the non-zero matrix el¬ 
ements, according to the second-order perturbation the¬ 
ory treatment of the second quantized Heisenberg hamil- 
tonian, between a localized and an itinerant electron spin 
coupled by Jf s . The term e(k + q) — e(k) corresponds to 
the energy diference between the ground state and the 
excited state. This function is known as a Lindhard func¬ 
tion that appears in the generalized form of the magnetic 
susceptibility 35 . Transforming the two summations into 
integrals, one may reach the actual expression for the 
RKKY interaction. Our calculated Fermi surface, shown 
in Fig. |9j evidences a rather complex dispersion relation 
for the GdCo 2 Zn 2 o compound, so obtaining an exact an¬ 
alytical solution from Eq. 10 is not possible. 

However, within a simple Fermi gas model (e(k) ~ k 2 ) 
the problem is simplified and more accessible. In this 
approximation the analytical expression of the RKKY 
interaction for the effective coupling between two lattice 
localized spins is given, in terms of k F , byE3 


Ag = 0.019(2) = J fs (0)r] Fs - 0.002(2), ( 8 ) 


b = 8 . 6 ( 6 ) Oe/K = —[F a (J%(q)) F r^ Fm ] + 1.9(6). (9) 
3dB 


Again from Eqs. 8 and 9 and using q Fs = 0.13(1) 
states/f.u., spin and eV, we obtain J/ s (0) = 167(7) meV 
and (J 2 s (g )) F 2 = 18(5) meV. Comparing this value with 
those for Jf p ( 0) and Jfd{ 0), we find that in these com¬ 
pounds the polarization component of the exchange pa¬ 
rameter Jf s ( 0 ) is an order of magnitude larger. 

Thus, in order to describe the antiferromagnetic or¬ 
dering of GdCo 2 Zn 2 o, it should be more appropriate to 
use an RKKY approach that considers only the exchange 
parameter J/ s (0) (due to the delocalized nature of the s- 
type electrons compared with p and d-type) rather than 
the Campbell mo del that considers the Jfd{ 0) term 

as the most important one. Therefore our obtained value 
in ESR analysis for J/ s (0) can be used to establish a cor¬ 
relation with the RKKY interaction, which is only valid 
in very dilute magnetic systems as this family of com¬ 
pounds. 

In general, the RKKY interaction depends strongly on 
the Fermi surface and can have different analytical forms 
for each case. Finding a suitable expression for a real 


^ 9ir(J /fl (0))W(r) 
eff ~ 64K F (2/c F ) 4 ’ 

( 11 ) 

where v is the number of conduction electrons per 
atom, Ep the Fermi energy, Jf s { 0) the coupling constant 
between the Gd 3+ spin Sf and the conduction electron 
spin s ce , {H = —Jf s (0)Sf . See), and F(r) is the RKKY 
function given by: 

sin( 2 fc F r) - 2k F r cos(2k F r) 

F{r) = - - A -. ( 12 ) 

Hence, in the GdCo 2 Zn 2 o compound, it is possible to 
estimate the effective exchange parameter, J e ff (Eq. 11), 
due to the coupling between two nearest Gd 3+ neighbors 
via the damped spin polarization of only the s-type ce, as 
it can observed in the third branch of the Fermi surface 
(Fig. [9] (c,g)). 

For GdCo 2 Zn 2 o the following band structure param¬ 
eters were determined from the DFT calculations: E F 
= 3.77 eV with k F = 0.99478 x 10 10 nr 1 [k F = 
y i \/2m e E F ) and v = 1.23 c.e/atom. Notice that the 
value of k F = 0.99478 x 10 10 m _1 is comparable to 
those reported for uncorrelated simple metals such as 
Cu, Ag and AuPS Then, the amplitude of the RKKY 
function F{r) at r = 6.0 A (nearest Gd-Gd neighbors) 
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gives F(6.0 A) = —7.99 x 10 37 m -4 (Fig. 11) and, conse¬ 
quently, J e ff = —25.2(7) x 10 -4 meV, a negative value 
as expected for antiferromagnetic ordering. 



r(A) 


FIG. 11: RKKY oscillatory coupling function F(r) as a func¬ 
tion of distance r between Gd ions. The arrow indicates the 
Gd-Gd separation in the GdCo2Zn2o system. 

Now, a microscopic Curie-Weiss temperature estima¬ 
tion can be made from this result: 37 


Oc 


2ZJ e ffS(S + 1) 
3 ks 


(13) 


with Z = 4 (Gd nearest neighbors in GdCo 2 Zn 2 o) and 
S = 1/2 for Gd 3+ . We obtain 9c = —1.2(2) K which 
is, within the accuracy of our experiments, in very good 
agreement with the bulk estimation (Fig. [7| extracted 
from the magnetic susceptibility measurements. 


(GdCo 2 Zn 2 o). The combination of ESR results with 
those of heat capacity, magnetic susceptibility and band 
structure calculations allowed us to estimate the polariza¬ 
tion component of the exchange parameters, Jf i(0 ) (i = 
s,p,d). Besides, by the assumption that only the s-type 
of ce can experience the bottleneck effect, due to their 
relatively weak spin-orbit coupling, we found that the 
average over the Fermi surface of the exchange parame¬ 
ter, associated to the ce momentum transfer, is different 
from the exchange parameter leading to local polariza¬ 
tion effects, i.e., (J/ s (9))f A Jfs{0 ). 

The exchange parameters obtained with this multi¬ 
band scenario revealed that J/ s (0) is dominant over 
Jf p ( 0) and Jfd{ 0). This allowed a tractable RKKY 
description for the antiferromagnetic behavior of the 
GdCo 2 Zn 2 o compound. Despite the fact that these com¬ 
pounds are structurally complex, we found that under 
certain reasonable approximations and using the combi¬ 
nation of different experimental results with DFT calcu¬ 
lations, the RKKY approach gave a very good quantita¬ 
tive description of the magnetic interaction, as expected 
for a naturally diluted structure of rare-earth ions. This 
was confirmed by the reasonably accurate prediction of 
the Curie-Weiss temperature in terms of microscopic pa¬ 
rameters. 

With the resulting establishment of GdCo 2 Zn 2 o as a 
model RKKY system, we expect that this work can 
provide key reference elements to help understand the 
behaviors of related materials such as RFe 2 Zn 2 o and 
YbT 2 Zn 2 o, with their more complex and remarkable elec¬ 
tronic and magnetic properties. 
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